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• Development of regional satellite-based pC02 algorithms for the North American east 
coast continental shelf 

• Assessment of the seasonal and interannual variability of surface ocean pC0 2 and sea-air 
CO 2 fluxes 

• Interannual estimates of the sea-air CO 2 flux show that the North American east coast 
continental shelf is a sink of atmospheric CO 2 ranging between 3.4 and 5.4 Tg C yr' 1 . 
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Abstract 


Underway and in situ observations of surface ocean pCOi, combined with satellite 
data, were used to develop pCOi regional algorithms to analyze the seasonal and interannual 
variability of surface ocean pCOi and sea-air CO 2 flux for five physically and biologically 
distinct regions of the eastern North American continental shelf: the South Atlantic Bight (SAB), 
the Mid- Atlantic Bight (MAB), the Gulf of Maine (GoM), Nantucket Shoals and Georges Bank 
(NS+GB), and the Scotian Shelf (SS). Temperature and dissolved inorganic carbon variability 
are the most influential factors driving the seasonality of pCOi- Estimates of the sea-air CO 2 flux 
were derived from the available pCOi data, as well as from the pCOj reconstructed by the 
algorithm. Two different gas exchange parameterizations were used. The SS, GB+NS, MAB, 
and SAB regions are net sinks of atmospheric CO 2 while the GoM is a weak source. The 
estimates vary depending on the use of surface ocean /?COi from the data or algorithm, as well as 
with the use of the two different gas exchange parameterizations. Most of the regional estimates 
are in general agreement with previous studies when the range of uncertainty and interannual 
variability are taken into account. According to the algorithm, the average annual uptake of 
atmospheric CO 2 by eastern North American continental shelf waters is found to be between 3.4 
and 5.4 Tg C yr' 1 (areal average of 0.7 to 1.0 mol CO 2 m' 2 yr' 1 ) over the period 2003-2010. 
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1.0 Introduction 


Coastal oceans, despite covering a small fraction of the earth’s surface, are important in 
the global carbon cycle because rates of carbon fixation, remineralization, and burial are much 
higher than the global average. A crucial difference between the coastal ocean and the open 
ocean is the proximity of sediments to the sea surface, providing a close coupling in space and 
time of the pelagic and benthic environments. Thus the shallow water column in coastal regions 
constitutes a close link between surface sediments and the atmosphere allowing relatively direct 
interactions between both the sedimentary and atmospheric compartments [ Borges et al., 2005; 
Thomas and Borges, 2012; Thomas et al., 2009; Thomas, 2004]. An additional characteristic of 
the coastal seas and continental shelves is the high temporal and spatial variability of CO 2 fluxes 
[Borges et al., 2005; Borges et al., 2008; Cai et al. , 2006; Frankignoulle and Borges, 2001; 
Shadwick et al., 2010; Shadwick et al., 2011]. The driving factors often vary within the system at 
seasonal time scales, and the deduction of general patterns remains difficult, typically requiring 
detailed case studies. 

The work of Borges [2005] was the first to compile a global coastal shelf sea-air CO 2 flux 
based on limited observed systems and using an up-scaling scheme. Borges [2005] showed that 
the inclusion of the coastal ocean increases the estimates of CO 2 uptake by the global ocean by 
57% for high latitude areas, and by 15% for temperate latitude areas, while at subtropical and 
tropical latitudes the contribution from the coastal ocean increases the CO 2 emission to the 
atmosphere from the global ocean by 13%. Cai et al. [2006] conducted a study of sea-air carbon 
exchange in ocean margins by grouping the numerous heterogeneous shelves into seven distinct 
provinces. Their results showed that the continental shelves are a sink of atmospheric CO 2 at 
mid-high latitudes (-0.33 Pg C yr' 1 ) and a source of CO 2 at low latitudes (0.1 1 Pg C yr' 1 ), with a 
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net uptake of -0.22 Pg C yr' 1 . Laruelle et al. [2010] evaluated the exchange of CO 2 between the 
atmosphere and the global coastal ocean from a compilation of sea-air CO 2 fluxes scaled using a 
spatially-explicit global typology of continental shelves. Their computed sink of atmospheric 
CO 2 over the continental shelf areas (-0.21±0.36 Pg C yr' 1 ) is at the low end of the range of 
previous estimates (-0.22 to -1.00 Pg C yr' 1 ). Laruelle et al. [2010] also concluded that the sea- 
air CO 2 flux per surface area over continental shelves, -0.7±1.2 mol CO 2 m' yr' , is twice the 
value of the open ocean based on the most recent CO 2 climatology at the time. More recently 
[ Cai , 2011] showed that the continental shelves are sinks of atmospheric CO 2 (-0.25 Pg C yr' 1 , 
but still with large uncertainty), accounting for -17% of open ocean CO 2 uptake (1.5 Pg C yr' 1 , 
Takahashi et al., 2009). The largest uncertainty of these scaling approaches stems from the 
availability of CO 2 data to describe the spatial variability, as well as to capture the relevant scales 
of temporal variability. 

Given that relatively large amounts of carbon are exchanged via the sea-air interface in 
coastal seas and continental shelves, the knowledge of the seasonal and interannual variability of 
the sea-air CO 2 flux in coastal oceans is a very important component of the carbon budget, which 
requires comprehensive regional studies. In general, the coastal ocean is characterized by a high 
variability in carbon cycling, which presents significant challenges in determining spatial and 
temporal integrals of relevant quantities, such as the sea-air CO 2 flux. Therefore, innovative 
methods are needed for scaling up relatively sparse field measurements, in this case surface 
ocean pCQj, into the required temporal and spatial resolutions to effectively derive regional sea- 
air CO 2 flux estimates. One method for obtaining such regionally integrated fluxes is through the 
use of biogeochemical-circulation models, which can be evaluated using the sparse field 
measurements, and then used to compute the mean and variability associated with these regional 
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fluxes [Hofmann et al, 201 1]. Satellite data, because of their high temporal and spatial 
resolution, provide another very promising asset to accomplish this goal. For example, Lohrenz 
and Cai [2006] conducted a satellite ocean color assessment of sea-air fluxes of CO 2 in the 
northern Gulf of Mexico. They used principal component analysis and multiple-regression to 
relate the surface ocean pC02 to SST, salinity and chlorophyll and used retrieval of 
corresponding MODIS-Aqua products to assess the regional distributions of pC02. 

In this paper we use multiple regression analysis to relate surface ocean pC02 to 
environmental variables (SST, surface salinity, and chlorophyll) and use the resulting equations 
with inputs from corresponding satellite products to provide an assessment of the spatial and 
temporal variability of the surface ocean /?COi and sea-air CO 2 flux for the North American east 
coast. A brief description of the biological/physical setting of the study region is provided in 
Section 2.0. The processing of in situ and satellite data sets and the development of regionally 
specific empirical pCC >2 algorithms are described in Section 3.0. The algorithm evaluation and 
the estimates of sea-air flux from the available pCC >2 binned data and algorithm are provided in 
Section 4.0, as well as a sensitivity analysis of parameters that influence the surface ocean / 2 CO 2 
seasonal and interannual variability. Finally, we provide a summary and discussion of suggested 
future work in Section 5.0. 

2.0 Physical and Biological Setting 

The temporal and spatial variability of the surface ocean pCC >2 on continental shelves are 
influenced by a combination of physical and biogeochemical factors, including surface 
temperature-driven solubility, biological processes, fall-to-winter vertical mixing, ocean 
circulation, river runoff, and shelf-ocean exchange [Wang et al., 2013], Here we provide a 
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summary of the physical and biological factors that are potentially important in shaping the pCOi 
variability in the North American east coast continental shelf. 

The definition of the coastal ocean is elusive, as it can be related to bathymetry, 
hydrography, or distance from shore; and some features, such as river plumes and coastal 
biomass maxima, can be ephemeral. Community efforts to standardize this definition to a fixed 
distance from shore, such as Hales et al. [2008] as adopted by the Surface Ocean CO 2 Atlas 
(SOCAT; http://www.socat.info/), extend seaward from the North American continent beyond 
what we feel represents the reach of coastal processes. As a result, we have used the outer 
boundaries of the regions defined by Hoffman et al. [2008, 201 1] to define the extent of the 
coastal ocean. The North American east coast (Fig. 1) encompasses three large regions of diverse 
physical and biological characteristics: the southeast U.S. continental shelf, also kn own as the 
South Atlantic Bight (SAB), the northeast U.S. continental shelf, and the Scotian Shelf (SS). 
Within the northeast U.S. continental shelf there are four sub-regions: the Middle Atlantic Bight 
(MAB), Georges Bank (GB), Nantucket Shoals (NS), and the Gulf of Maine (GoM). For this 
study we combined the GB and NS regions into a single region (GB+NS) for simplicity and 
because these two regions share many similar physical and biogeochemical attributes [Fox et al., 
2005; Shearman and Lentz, 2004; Thomas et al, 2003]. These North American continental shelf 
sub-regions are defined in Hofmann et al. [2011] with the GB+NS region separated from the 
GoM as in Hofmann et al. [2008]. The 58 coastal sub-regions shown in Hofmann et al. [2008] 
were developed based on a combination of bathymetry, SST fronts, stratification, and biological 
properties. For simplicity, here we consolidate the very fine regional domains into five major 
sub-regions described above. However, we recognize that previous studies have adopted other 
methods to identify regional domains [ Hales et al., 2008; Hales et al., 2012]. For example, a self 
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organizing mapping method has been adopted to sub-regionalize the North American Pacific 
Coast [Hales et al., 2012]. The method relies on an artificial neural network to identify 
biogeochemical regions within the target study area. 

Our focus is on the continental shelf which we operationally define as depths less than 
200 meters since the depth of the actual shelf break varies. Bathymetric variation in our study 
area is large. Portions of GB and NS are only several meters below the sea surface, whereas in 
the GoM and areas of the SS, water depths exceed 200 m. Our study area is also at the 
‘crossroads’ of the north-flowing Gulf Stream and the southwest-flowing slope water-Labrador 
current [Rossby, 1987]. Chapman and Beardsley [1989] suggest that glacial melt and runoff 
from Western Greenland generates a buoyancy-driven coastal current that flows over the SS and 
GB and eventually into the MAB. This coastal current is an important driver to the distribution 
of the marine CO 2 system, including surface pCOi along its flow path [Wang et al., 2013], i.e., 
the Gulf of St Laurence, the SS, the GoM and the MAB. There is little exchange of water 
between the MAB and SAB along the narrow shelf at Cape Hatteras. In the SAB, the Gulf 
Stream is close to the shelf break and has a direct influence on the outer SAB shelf [Signorini 
and McClain, 2007], readily identifiable by the warm and salty signature shown in seasonal 
maps of sea surface temperature (SST), sea surface salinity (SSS), and chlorophyll (Chi) of Fig. 
2 (see Section 3.0 for methodology), whereas north of Cape Hatteras, the influence of the Gulf 
Stream is more indirect. Here anti-cyclonic warm core rings result from landward meanders of 
the Gulf Stream [Joyce et al. , 1992]. The rings are carried in the southwestward flow of slope 
water where they interact with the outer shelf from GB to Cape Hatteras, frequently entraining 
phytoplankton-rich shelf water [Joyce et al, 1992]. Near Cape Hatteras, the warm core rings 
may be reabsorbed into the Gulf Stream, a process readily apparent in daily time series 
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animations of chlorophyll (Chi) and SST. In the SAB, the outer shelf waters are warmer (Fig. 2) 
in summer and autumn than winter and spring due, in part, to the proximity of the Gulf Stream as 
a result of the expansion of the subtropical gyre [ Signorini and McClain, 2007]. 

The pCOi variability in riverine-plume systems is a result of complex biogeochemical 
interactions. In the Gulf of Maine for instance, labile riverine carbon is responsible for sustaining 
supersaturated pCC >2 conditions in late fall, while at other times of the year phytoplankton 
productivity, most likely driven by inputs of riverine dissolved inorganic nitrogen, is responsible 
for pCC >2 undersaturation [Salisbury et al., 2008]. The North American east coast continental 
shelf is influenced by the discharge of several major rivers and estuaries (Chesapeake Bay, 
Delaware Bay, and Gulf of St Lawrence, for example) that contribute to complex physical and 
biogeochemical interactions that influence the seasonal and interannual variability of the surface 
ocean pCC> 2 , an important parameter for the determination of the sea-air CO 2 flux. Vandemark et 
al. [2011] showed that the observed pC02 and CO 2 flux dynamics in the Gulf of Maine are 
dominated by a seasonal cycle, with a large spring influx of CO 2 and fall-to-winter efflux back to 
the atmosphere. They also showed that in the western Gulf of Maine the ocean is a net source of 
carbon to the atmosphere (+0.38 mol CO 2 m" yr" ) over a period of five years, but with a 
moderate interannual variation where years 2005 and 2007 represent cases of regional source 
(+0.71 mol CO 2 m" yr' ) and sink (-0.1 1 mol CO 2 m" yr' ) anomalies, respectively. Comparison 
of results with the neighboring Middle Atlantic and South Atlantic Bight shelf systems showed 
that the Gulf of Maine differs by enhanced pC02 control factors other than temperature-driven 
solubility, such as biological drawdown, fall-to-winter vertical mixing, and river runoff 
[Salisbury et al., 2008; Shadwick et al, 2010]. 
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Shadwick et al. [2011] investigated the seasonal variability of pCOo in the Scotian Shelf 
and concluded that the region acts as a net source of CO 2 to the atmosphere on an annual basis 
(1.4 mol CO 2 m' yr" ). On a seasonal basis, there is a reversal of the flux only when a 
pronounced undersaturation of surface waters is reached for a short period during the spring 
bloom. Outside of the spring bloom period, the competing effects of temperature and biology 
influence on surface pCOi are nearly equal and opposite. DeGrandpre et al. [2002], based on 
measurements of surface ocean pCOi during the Ocean Margins Program [ Verity et al, 2002], 
concluded that the MAB is a sink of atmospheric CO 2 with an annual mean of -1.0 ± 0.6 Tg C yr' 
, or an area average of -1.1 ± 0.7 mol CO 2 m' yr' . A significant portion of this atmospheric 
uptake is a result of the annual cycle of heating and cooling combined with strong winds during 
the winter undersaturation period. 

Jiang et al. [2008] showed that on an annual basis the SAB is a relatively small net si nk 
of atmospheric CO 2 (-0.48 ± 0.21 mol CO 2 m' 1 yr' 1 ). Seasonally, the SAB shifts from a sink of 
atmospheric CO 2 in winter to a source in summer. The annual cycle of sea surface temperature 
plays an important role in controlling the seasonal variation of pCOi- The combination of 
stronger wind speeds during fall-winter, when CO 2 undersaturation is significant due to lower 
SSTs, results in a net annual CO 2 sink. Other important factors controlling the pCOi variability 
in the SAB are the marsh export of organic carbon and DIG in the wann months (June- 
November), which directly supports CO 2 outgassing in these months via organic carbon 
decomposition and increase in DIC [Jiang et al., 2013; Wang et al., 2005]. I 11 addition, the marsh 
areas in the SAB also export alkalinity, another important factor influencing the variability of 
/ 7 CO 2 and sea-air flux [ Wang et al., 2005; Wang and Cai, 2004]. 
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The seasonal Chi climatology from MODIS Aqua (Fig. 2) shows that the maximum Chi 
in the GoM, GB and NS occurs during spring (March-April-May, MAM). The GB region has the 
highest Chi in spring, but it is maintained at concentrations above 2.5 mg m' in all seasons due 
to vigorous tidal mixing. Fig. 2 also shows that the low-salinity nearshore waters along the entire 
east coast coincide with regions of elevated Chi, an indication of the influence of nutrient-rich 
riverine waters. On the MAB shelf, there is a high-C/ 7 / region during winter (December-January- 
February, DJF) in the near shore and outer-shelf waters, but the fall bloom (SON) dominates 
between approximately the 40- and 60-m isobaths. The high satellite-derived ‘C/ 7 /’ in winter may 
be in part colored dissolved organic matter flowing out from rivers, plus photo-acclimation by 
phytoplankton (higher Chl-a due to low surface solar radiation and a well-mixed water column). 

The minimum surface Chi over much of the MAB occurs during summer (JJA) when 
highest SST (Fig. 2), peak stratification and a pronounced subsurface Chi maximum layer occur 
[O'Reilly and Zetlin, 1998]. Summer mixed-layer depths of ~3.5 to 10 m are typical for MAB 
shelf waters. The spring bloom (MAM) is clearly shown by the elevated Chi concentrations in 
the MAB, GB, and GoM (Fig. 2). Fig. 2 also shows that the SAB Chi has its largest changes in 
the outer shelf, with a maximum in DJF and lowest values in JJA under the influence of the 
oligotrophic waters of the Gulf Stream. 

3.0 Data Sets and Methods 

3.1 Processing of In Situ and Satellite Data Sets 

The surface ocean pCOi data are obtained from SOCAT, combined with additional 
available data from regionally specific field experiments (see Appendix A) and binned by month 
for each year (1978-2010) into 0.15°x0.15° grid cells. The SOCAT data [Pfeill et al, 2012] holds 
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6.3 million quality-controlled surface ocean pCOi from the global oceans and coastal seas 
covering the period of 1968 to 2007. These data were put together following unifonn format and 
a strict protocol that included quality control with clearly defined criteria performed by a team of 
international experts. 

The MatLab function bin2d, developed by J. Nielsen and available at the Nansen 
Environmental and Remote Sensing Center (NERSC) from www-2.nersc.no/~even/, was used to 
bin all data sets into the study grid. First, all the available data within 24°N to 46°N and 82°W to 
56°W were selected for binning. These included 416,261 co-located surface ocean /?CC> 2 , SST 
and sea surface salinity (SSS) values from SOCAT from the period 1978-2007, 1 1,628 from the 
2006 SAB cruise (only 2005 cruises are included in SOCAT), and 309,665 from the GoM (2004- 
2010). The binned pCOo data were then adjusted to reference year 2004 using an atmospheric 
growth rate of 1.68 palm yr' 1 [Le Quere et al., 2010] and assuming that the surface ocean pCOi 
is trending at the same pace as the atmosphere. All the adjusted pCXX data were then binned into 
12 individual calendar months, each containing the average of all data within a particular month 
and grid bin. The data were then divided into regional study domains following the boundaries 
shown in Fig. 1. 

The available pC02 data were divided into two individual sets, one dedicated to 
algorithm development (data bins covering more than 6 months) and one dedicated to algorithm 
evaluation (data bins covering less than 6 months). Surface ocean pCOi data from underway 
(UW) transects across the Scotian Shelf and pCOj time series from the CARIOCA buoy located 
at 44.296°N and 63.257°W [ Shadwick et al., 2010] were also used for algorithm evaluation, 
together with SOCAT data on the Scotian Shelf not used for the algorithm development. Fig. 3a 
shows color-coded SOCAT surface ocean /?COi cruise tracks and Fig. 3b shows corresponding 
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coastal binned data with associated color-coded temporal coverage in months. The highest 
temporal coverage corresponds to the most travelled routes (in orange to red), i.e., most frequent 
destination ports (Boston, New York, Norfolk, Miami) used by the Volunteering Observing 
Ships (VOSs). The VOS ships according to map available at the CDIAC web site 
(http://cdiac.ornl.gov/oceans/VOS_Program/) are: the Skogafoss, A. Companion, Oleander, 
Falstaff, and Explorer of the Seas. The SOCAT data set also includes transects occupied by 
research vessels. Fig. 3 clearly shows that the surface ocean pCOi data have spatial and temporal 
distribution gaps that may be potentially responsible for biases in the calculation of sea-air 
fluxes. 

Monthly sea-surface salinity (SSS) climatology was interpolated and gridded onto the 
0.15°x0.15° study domain grid using the World Ocean Database (WOD) 2009 station data and 
the method of Kriging. The Interactive Data Language (IDL) function KRIG2D was used for this 
purpose. Monthly climatologic mixed layer depth (MLD) was derived from WOD 2005 for the 
entire East Coast based on temperature profiles using 0.5° C temperature difference criterion 
[ Hofmann et al, 2008]. The MLD data were binned into the same 0.15°x0.15° study domain 
grid. 

Both data and algorithm sea-air CO 2 flux estimates were obtained using gridded 
(0.25°x0.25°) winds from the Jet Propulsion Laboratory Cross-Calibrated Multiple Platforms 
(CCMP, Atlas et al., 2011) product (ftp://podaac-ftp.ipl.nasa.gOv/allData/ccmp/L2.5/flk) . 
Monthly wind climatology was derived using data from 1999 to 2008, a period approximately 
centered on 2004, the reference year adopted for the adjusted surface ocean pCO^ data. The 
climatologic and interannual CCMP monthly winds were re-gridded (0.15°x0.15°) and 
extrapolated nearshore using the function “surface” from Generic Mapping Tools (GMT, Smith 
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and Wessel, 1990; Wessel and Smith, 1991) which is based on an adjustable tension continuous 
curvature surface gridding method. High frequency (10-minute) winds from 10 NOAA National 
Oceanographic Data Center NDBC buoys (http://www.nodc.noaa.gov/BUOY/) and hourly winds 
from Sable Island were used to obtain correction coefficients to account for nonlinearities in the 
gas exchange parameterization resulting from the use of monthly climatologic winds. The 
method for deriving these coefficients is described under sub-section 3.3. 

All parameters used to develop the pCOi algorithm and to derive the sea-air CO 2 flux, 
including all satellite data products (SST and Chi), SSS and the CCMP wind speed were also 
binned monthly into the same grid. The satellite data products consisted of 9-km, level 3 
mapped, MODIS Aqua (MODISA) climatologic and interannual monthly composites of SST and 
Chi obtained from the NASA ocean color distribution archive (http://oceancolor.gsfc.nasa.gov/) . 
A validation between log-transformed MODISA Chi retrievals vs. all available in situ 
observations (SAB to GoM, depth<=200m, ,V=404), conducted using the SeaBASS (SeaWiFS 
Bio-optical Archive and Storage System: http://seabass.gsfc.nasa.gov/) data search and 
validation tools, showed good matchup agreement (r 2 =0.75, RMSE=0.30, APD=35.8%). For the 
algorithm development we used the available binned surface ocean pC02, SST and SSS derived 
from the in situ data, combined with monthly climatologic satellite Chi binned at the same grid 
points as no in situ concurrent Chi measurements are available. For the algorithm application we 
used monthly interannual (2003-2010) satellite SST and Chi, and monthly climatologic SSS 
derived from WOD 2005 data. 

Seasonal maps were constructed by averaging the monthly data and derived products into 
four three-month composites, defined as: winter (December-January-February, DJF), spring 
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(March-April-May, MAM), summer (June -July-August, JJA), and autumn (September-October- 
November, SON). 

3.2 Development of Regional pC0>2 Algorithms 

The algorithm development is based on binned in situ pCOi, SST and SSS, and satellite- 
derived Chi monthly climatology, as well as day of the year (Julian day). The algorithm was 
developed through the multiple linear regression (MLR) analysis based on all spatial bins 
containing more than six available monthly occurrences of the in situ data (remaining data were 
reserved for evaluation), and is represented as: 

pC0 2 — [a + bDay' + c(T — T 0 ) + d{S — S 0 ) + e[log 10 (Chl ) — log 10 {Chl o )]+ 

\.68{year - 2004) (1) 

/2n{Day - y)\ 
where, Day = cos ( — J 

The first terms in brackets represent the surface ocean /?COi corrected to the year 2004 and the 
last term is a correction factor for different years to account for the rise of surface ocean p CO 2 
due to the uptake of anthropogenic CCL. The input for “ Day ” (Julian day) was normalized 
sinusoidally ( Day ') to emphasize the seasonal cycle and to allow January to be close to both 
February and December [Friedrich and Oschlies, 2009; Lefevre et al., 2005]. The value of y 
(phase of Day" in days) is optimized via iteration (ranging from 0 to 365 days) until the 
minimum RMSE is obtained. T 0 , S 0 , Chl Q are temperature, salinity, and chlorophyll mean values 
for each region. The choice of logio {Chi) instead of Chi in our algorithm was an arbitrary choice, 
and therefore limited mechanistic infonnation can be drawn in the empirical result. 
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A separate analysis was conducted to evaluate the algorithm by using surface ocean pC0 2 
data not used in the development of the algorithm equations (see Section 4.1). These data 
consisted of bins from the monthly composites that have less than six months of available pC0 2 
occurrences. Satellite-derived SST, Chi, in situ SSS monthly climatology was matched with the 
locations and months of the selected pCOi bins and used as algorithm input. The pC0 2 derived 
from the algorithm (pCO^) was matched with the observed pC0 2 ( pCO^ ) and a scatter plot 
and histogram of residuals were made for all combined regions to evaluate the algorithm 
perfonnance. The algorithm was also evaluated using data from the SS ( Shadwick et al., 2010). 

3.3 Calculation of the Sea-air CO 2 Flux 

The sea-air pC0 2 difference (ApCCT) was calculated using monthly GLOBALVIEW 
[GL0BALVIEW-C02, 2011] atmospheric xCCV from Grifton, North Carolina, a station located 
approximately midway in the study domain. The xCCV (in pmol mol’ 1 ) was converted to pC0 2 
(air) using the method of Jiang et al. [2008]. For this conversion we used monthly surface 
barometric pressure and air temperature from NOAA NCEP-NACR CDAS-1 [Kalnay et al., 
1996] and monthly climatologic SSS from WOA09. Although several other GLOBALVIEW 
stations are available along the study coastal domain, the atmospheric pC0 2 records are not very 
different to justify a more site-specific use of the data. Regarding the use of the atmospheric 
XCO 2 in this study, it has been demonstrated that there are uncertainties involved in using marine 
boundary layer XCO 2 rather than the in situ xCCV due to the effect of continental processes. For 
example, Jiang et al. [2007] showed that the average atmospheric XCO 2 on the SAB can be 
almost 10 ppm higher than the measured in the open ocean with the potential of reversing the 
direction of the sea-air flux. Although this is a potential source of uncertainty in the calculation 
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of the sea-air flux, concurrent in situ atmospheric xCOi are only available for a limited number 
of coastal cruises. 


Climatologic (1999-2008) CCMP monthly wind speeds at 10-m anemometer height 
(£/io), based on a decade of data centered on the reference year 2004, were binned similarly and 
used to derive the monthly sea-air CO 2 flux for each bin and each month using the following gas 
transfer parameterization 


Flux — /c 660 (~) -1 ^ 2 s ApCO z (2) 

2 1 

in units of mol CO 2 m" d' . Sc is the Schmidt number (non dimensional), s the solubility of CO 2 
in seawater in mol CO 2 m' palm" , and A pC0 2 is the sea-air pC0 2 difference in palm. The term 
&660 is the quadratic gas transfer coefficient in cm If 1 (converted to m d' 1 ). We calculated the sea- 
air CO 2 flux using two relationships of gas exchange with wind speed (C/ 10 ), the quadratic 
dependence formulation of Ho et al. [2011], for which k 660 — 0.2 62 C 2 t/ 1 2 0 , and the polynomial 
dependence of Wanninkhof et al. [2009], for which k 660 — 3 + 0.1 U 10 + 0.064C 2 t/io + 
0.011C ?1 U 20 , using the appropriate nonlinearity correction coefficients C 2 and C 3 , which are 
correction factors to account for the use of monthly climatologic wind speeds [ Jiang et al, 

2008]. These were calculated using 10-minute wind speeds from 10 NDBC buoys distributed 
within the SAB, MAB, GB+NS, and GoM regions, and Sable Island 1-hour wind speeds for the 
SS (see locations in Fig. 1), and the correction factor equations given in Jiang et al. [2008] , 

C -2 = Qz" Uj) / U mean and C 3 = Qs” I//) l^mean, where Uj is the high-frequency wind 
speed (m/s), U mean is the monthly mean wind speed (m s" 1 ), and n is the number of available wind 
speeds in each month. The value of C 2 and C? were obtained for each site and month for the 
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period 1999-2008. Monthly climatologic averages were calculated for each site and for each 
region. The values of C 2 range from 1.2 to 1.3, while those for C 3 range from 1.6 to 2.0. These 
values were then used to apply corrections to the gas transfer parameterizations when calculating 
the sea-air CO 2 flux. The same methodology was applied to derive data-based and algorithm- 
based sea-air fluxes. We use the atmospheric convention for the CO 2 flux, i.e., a negative flux is 
defined as a sink of atmospheric CO 2 by the ocean. 

The regional algorithms (Table 1 and Equation 1) were used to derive values of surface 
ocean pCC >2 using MODIS Aqua monthly composites of SST and Chi for 2003-2010, and 
monthly SSS climatology. Gap filling of missing satellite data was done with monthly 
climatology composites for each of the input parameters. The sea-air CO 2 flux was then 
computed using interannual monthly CCMP winds and the gas transfer parameterization shown 
in Equation 2. 

3.4 Monthly Climatology of DIC and Alkalinity for pCOi Parameter Sensitivity 

The data sets used to generate monthly climatologies of DIC and alkalinity (Aik) include 
the MODIS SST monthly climatology, the Kriged monthly SSS climatology derived from WOA 
2009 salinity data, and surface ocean pCOi from the algorithm. Monthly alkalinity was derived 
as a function of salinity from Cai et al. [2010] using SSS monthly climatology. DIC was then 
derived from alkalinity, SST, SSS, and monthly pCOi from the algorithm using C02SYS 
(http://cdiac.oml.gov/ftp/co2sys/C02SYS calc MATLAB/) , a MatLab program to calculate the 
state of the carbonate system. The input for C02SYS consisted of alkalinity, DIC, SST, SSS, the 

choice of H 2 C0 3 and HCO s dissociation constants (AT 1 , Ki) of “Mehrbach refit” [Dickson and 

Millero, 1987], the choice of HS0 4 dissociation constant of “ Dickson” [Dickson, 1990], and 
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zero concentration for silicate and phosphate. The total borate-salinity relationship of Uppstrom 
[1974] was used. 

The monthly binned SST, SSS, DIC, and alkalinity fields were then averaged over each 
region to obtain 12 monthly values for each variable and region. We refer to these regional 
averages as SST 1 , SSS 1 , DIC 1 , and Aik 1 , where the superscript indicates the calendar month from 
1 to 12. We also computed the annual average of each of these four spatial averages, which we 
call, SST, SSS, DIC, and Aik. From the regional averages, we computed the monthly pCOi using 
C02SYS, 

pCOi, = pC0 2 (SST 1 , SSS j , DIC’, Aik’), (3) 

and the annual average, pC0 2 . 

The deviation of pCO 2 from its annual average is given by 

8 l = pCO], - pC0 2 (4) 

To detennine the sensitivity of pCXh to each of the four variables, we hold three variables 
at their annual averages and let the fourth variable change from month to month. For example, to 
detennine the impact of temperature on pC0 2 , we computed 

pCOf T = pC0 2 (SST i ,SSS7DIC, Aik) (5) 

In an analogous way, we computed pCCC' 855 , pCCF 1 ’ 1 ^ , and pCCF'" 811 ) which describe 
the respective influences of SSS, DIC, and Aik on /?C0 2 . We also computed the deviation of 
pCCF from its annual average due to each of the four variables. For example, the deviation of 

pCC >2 from its annual average due to temperature is 5 1, = pCO^ — pC0 2 . Similarly, 8 ’ , 
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S 1,DIC , and S 1,Alk , describe the deviations of pC0 2 from its annual average due, respectively, to 
SSS, DIC, and Aik. The results of this analysis will be discussed in sub-section 4.3. 


4.0 Results and Discussion 

Regional algorithms were developed with distinct coefficients derived for each of the five 
regions (Table 1) and then used to derive seasonal and interannual surface ocean y>COi and sea- 
air CO 2 fluxes (Tables 2 and 3). 

4.1 Perfonnance of Regional Algorithms 

In this section we provide an assessment of the statistical importance of each proxy 
parameter used in the algorithm (Fig. 4), regional matchups of algorithm versus data and 
seasonal pCOi plots based on monthly averages derived from data and algorithm (Fig. 5), 
algorithm versus data matchups using pC0 2 observations not used in the algorithm development 
(Fig. 6), a regional matchup analysis for the Scotian Shelf (SS) using a combination of UW pCOi 
data from Dalhousie University and a few from SOCAT (Fig. 7), and time-series of algorithm 
pC0 2 for seven distinct sub-regions of the SS (concurrent data points) following a more recent 
work of Thomas et al. [2012] (Fig. 8). Finally, a high frequency algorithm validation was 
perfonned against surface pCO >2 observations from the CARIOCA buoy on the SS using 
concurrent hourly observations of SST, SSS, and Chi (Fig. 9). 

Fig. 4 shows the statistical (goodness-of-fit) perfonnance resulting from the incremental 
addition of proxy parameters for each of the five regions. The statistical perfonnance is shown as 
a goodness-of-fit diagram with normalized RMSE on the x-axis, and (1 - r ) on they- ax is. 
Consequently, a perfect fit would lie at the origin of this diagram (0, 0). The diagram shows that 
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2 

the variable Day ’ by itself provides (1 - r) values less than 0.6 for all regions. Incremental 
improvements of both normalized RMSE and (1 — r ) are different for each region. Extreme 
examples of statistical improvement are the addition of salinity for the SAB and logio (Chi) for 
the SS. 

Fig. 5 shows scatter plots of algorithm-derived versus observed surface ocean pC0 2 and 
associated seasonal plots of regionally- averaged pC0 2 . As shown in Table 1, there is a statistical 
range for the coefficients derived for each region using Equation 1 . The r is lowest for the GoM 
(0.42) and highest for the SAB (0.82). The quality of the statistical fit depends on a combination 
of factors, including data coverage and how well the proxy variables represent the surface ocean 
pC0 2 variability in space and time within each region. 

The regional algorithms were then applied using binned inputs (SST, SSS, and Chi ) 
matching the month and location of the observed surface ocean pC 02 not used for the algorithm 
development, and then compared with the corresponding observed pC0 2 . The results are shown 
in Fig. 6a and 6b. The observed versus algorithm correlation coefficient (color coded scatter plot 

9 

in Fig. 6a with summary of statistics in the legend) range from 0.27 (r ) for the GoM with a 
RMSE= 25 patm to 0.78 for the SAB with a RMSE=21 patm. The histogram of residuals (Fig. 
6b) shows that 86% of the residuals are less than the observed pC0 2 standard deviation (±a), 
while 40% of residuals are within less than a/3 (±16 patm). 

Data from SOCAT on the SS, and Dalhousie University UW transects [ Shadwick et al, 
2010] covering the period of 2004-2008, were averaged within seven 2°x2° boxes on the SS (Fig. 
7a) and compared with area-averaged algorithm predictions within the same boxes. The scatter 
plot of observed vs. algorithm pC0 2 for the 37 resulting averages is shown in Fig. 7b. The 
agreement between data and algorithm predictions is quite reasonable with r =0.79 and 
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RMSE=26.2 |iiatm. The time series of algorithm pC0 2 was obtained using SST and Chi from 
MODIS Aqua monthly composites and WOA09-derived SSS climatology. The algorithm time 
series for all seven boxes are shown in Figs. 8a and 8b with the SOCAT (red circles) and UW 
(blue circles) values superposed for comparison. A high frequency algorithm test was done by 
comparing the CARIOCA buoy one-hour pC0 2 record on the SS during 2007-2010 with 
algorithm results using one-hour inputs of SST, SSS and calibrated fluorometer Chi concurrent 
observations from the buoy. These data have been reported by Thomas et al. [2012]. The time 
series and scatter plot of observed vs. algorithm pC0 2 are shown in Fig. 9. The algorithm 
predictions track the observed pC0 2 reasonably well with r =0.46, RMSE=40.3 patm and mean 
absolute percent difference (MAPD) of 8.8%. The observed and algorithm values for 2007-2010 
mean and standard deviation are quite similar, 422.3±54.7 palm and 413.1±56.9 patm, 
respectively, which show a relatively small bias (9 patm) and very similar variance. 

4.2 Seasonal Surface Ocean pC0 2 , Alkalinity, DIC and Sea-Air Flux from Data and Algorithm 
Fig. 10 shows seasonal maps of algorithm surface ocean pC0 2 adjusted for reference year 
2004 and corresponding seasonal maps of alkalinity and DIC. Fig. 10 shows that the temporal 
and spatial variability of pC0 2 is quite different from region to region and that the seasonal 
changes are not in sync among the five analyzed coastal domains. This is also evident in the 
seasonal plots of data-derived surface ocean pC0 2 in Fig. 5. The lowest values (280 to 320 patm) 
occur mostly during winter (DJF) in the MAB, SAB, and in the near shore areas of the SS in 
spring (MAM). Low values are also present in spring in the GB+NS region. These low values are 
generally associated with low SSTs (See Fig. 2). The highest values (> 480 patm) occur in the 
offshore region of the SS in autumn (SON) and the near shore areas of the SAB in summer 
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(JJA), the latter influenced by the discharge of carbon-rich (primarily DOC) estuarine effluents 
[Alberts and Takacs , 1999; Cai, 2011] and marsh DIC export [Wang and Cai, 2004]. The surface 
ocean pCOi in the MAB shows much less variability alongshore than cross-shelf, except in the 
southern region and outer shelf where Gulf Stream intrusions and shelf-slope fronts induce 
strong hydrographic and biogeochemical horizontal gradients. DeGrandpre [2002], and 
references within, identified similar alongshore homogeneity in connection with little alongshore 
variability on mid-shelf hydrography, nutrients, surface dissolved oxygen, Chi concentrations, 
and primary production. The high values in the offshore region of the SS in autumn are 
associated with low drawdown by phytoplankton, as indicated by the higher values of DIC, as 
shown in Fig. 10 discussed later in this section, and confirmed by the work of Craig et al. [2013] 
for this region. The GoM has highest pCC >2 (> 400 palm) values in winter and fall when vertical 
mixing is more vigorous and phytoplankton drawdown is significantly reduced. 

The seasonal maps of alkalinity in Fig. 10 follow the seasonal surface salinity distribution 
in Fig. 2 as alkalinity was derived as a linear function of salinity, albeit with different 
coefficients for each region. There is a sharp transition in alkalinity at Cape Hatteras. South of it, 
in the SAB, alkalinity is highest in the middle and outer shelves due to the influence of high- 
salinity Gulf Stream waters. Alkalinity is highly reduced in the nearshore region under the 
influence of low-salinity riverine waters. However, in the very nearshore areas high alkalinity 
values were observed due to significant export from the marsh areas during the warm months 
[Cai et al, 1998]. North of Cape Hatteras all regions have much lower alkalinity than the middle 
and outer shelf regions of the SAB. The inner and middle shelf regions of the MAB and southern 
GoM have even lower alkalinity, especially during summer (JJA) when surface salinity is at a 
minimum. This summer minimum salinity follows the peak discharge of the major rivers in 
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spring with a delay of approximately 1-2 months [ Whitney , 2010], However, the SSS minimum 
on the SS comes in autumn (SON) with the peak St Lawrence outflow. 

The Aik and salinity relationships generally followed a single river-ocean mixing line in 
the SAB and MAB regions, but a two-segment line in the northeastern waters due to the strong 
alongshore current and influences from the low alkalinity local rivers [Cai et al., 2010]. 

The seasonal DIC maps in Fig. 10 show highest values in the GoM and offshore regions 
of the SS in winter-spring, a likely result of vigorous vertical mixing. Lowest DIC values occur 
in the MAB and southern GoM in summer, influenced by the low-DIC riverine waters that peak 
during spring, as well as low-DIC water of the Labrador Coastal Current that flows through the 
region [Wang et al., 2013]. The DIC seasonal variability is also highly influenced by the 
drawdown of CO 2 by the net community production during spring-summer. In general, the SAB 
has much less seasonal DIC and alkalinity variability than the other regions to the north. 

The monthly and annual mean sea- air CO 2 flux was calculated using ApC0 2 derived from 
both binned data and algorithm (Table 2) and the two gas transfer parameterizations described in 
Section 3.3. The estimates were based on monthly wind climatology for 1999-2008 derived from 
satellite (Atlas CCMP) winds. The differences between the two different parameterizations are 
relatively small ranging from 6% to 17%, except for the GoM where the fluxes are small causing 
much larger differences between the two methos. For simplicity we compare the flux estimates 
between binned data and algorithm based on the Ho et ah [2011] parameterization. 

There is a general agreement in sign and magnitude between the data-derived and 
algorithm-derived estimates for the MAB, SAB, and GB+NS (Table 2). The annual mean sea-air 
CO 2 flux in the GoM derived by both methods range from +0.02±0.12 to +0.17±0.32 Tg C yr' 1 , 
or a weak source to the atmosphere on average, but within the range of the estimates given by 
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Vandermark et al. [2011] for the southern GoM (-0. 16 to +1.1 Tg C yr' 1 when converted from 
specific to up-scaled total sea-air flux for the entire GoM). The MAB, SAB, GB+NS and SS are 
net si nk s ranging from -0.6±0.2 to -1.8±0.2 mol CCfi m' yr' . These estimates from the binned 
data and algorithm are in general agreement with previous studies (see Table 2) when the range 
of uncertainty and interannual variability are taken into account. One exception is the SS where 
previous studies [Shadwick et al., 2010; Shadwick et al., 2011] indicate that the SS is a source of 
CO 2 to the atmosphere while this study indicates the opposite. Since the algorithm seems to 
perform well in the SS when compared with the available data, the reason(s) for the apparent 
discrepancy remains elusive and highlights the fact that there are still large differences in the sea- 
air flux estimates with different degrees of uncertainty from region to region. 

The combined uptake by the east coast continental shelf based on both binned data and 
algorithm, and using both gas transfer parameterizations, ranges from 3.6 to 4.3 Tg C yr' 1 . 

4.3 Sensitivity Analysis of Parameters that Influence the pCXh Seasonal Variability 

Here we present a sensitivity analysis of the most influential parameters affecting the 
surface ocean pCOi variability in the study region. The seasonal cycles of each influential 
parameter are plotted in Fig. 1 1 together with the seasonal surface ocean pCOj from the 
algorithm with the seasonal mean removed. Inspection of Fig. 1 1 shows that the amplitude of 
SST and DIC contributions in the MAB, GoM, GB+NS, and SS are similar but having opposite 
phase. Seasonal variability of pCOi (DIC) in these regions is consistent with winter mixing 
enhancement and biological drawdown in spring-summer. In contrast, the major contributing 
factor to the seasonal pC0 2 variability in the SAB is SST. Alkalinity influence is the third most 
important and salinity relatively the least influential. However, salinity has an impact in the 
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statistical improvement of the pCOi algorithm, most pronounced in the SAB, which is a region 
where seasonal SSS variability is large (see Fig. 2), especially on the inner shelf. 

The seasonal DIC variability averaged for all five study regions, with the MLD 
superimposed, is shown in Fig. 12. The four study regions north of Cape Hatteras (MAB, GoM, 
GB+NS, and SS) have distinct DIC seasonal cycles with amplitudes of 100 to 120 pmol kg' 1 . 
Regionally-averaged winter MLDs range from 30 m in the MAB to more than 100 m in the 
GoM. Deeper MLDs in winter/autumn, resulting from wind and convective mixing, is the major 
factor contributing to the elevated DIC concentrations (2010 to 2080 pmol kg' 1 ) shown during 
these seasons. The shoaling of the MLDs in spring-summer, together with the drawdown of CO 2 
by biology, are the major factors driving the significant reduction in surface DIC. For instance, in 
the MAB the DIC drops from 2020 pmol kg' 1 in February-March to 1900 pmol kg' 1 in June. In 
addition to biology and deep mixing, DIC, and consequently the surface ocean pCOi, is also 
affected by sea-air exchange. In the GoM, for instance, there is a significant effect of the sea-air 
exchange on DIC when the SpCOi is high and the mixed layer becomes very shallow (J. 
Salisbury personal communication, 2012). The amplitudes of the seasonal MLD and DIC in the 
SAB are significantly less than in the other regions, most probably due to the shallower depths 
and much lower phytoplankton productivity. 

4.4 Interannual Variability of Surface Ocean pCC >2 and Sea- Air Flux 

The interannual variability of surface ocean pCC >2 and sea-air CO 2 flux were calculated 
using the algorithm (Equation 1) with inputs from monthly satellite products (SST and Chi) for 
2003-2010 and climatologic SSS. The sea air flux was computed using monthly CCMP winds 
for the same period. The results are shown in Fig. 13 (pC02 left panel, sea-air flux right panel) 
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and summarized in Table 3. Note that the algorithm results in Table 2 were derived using 
monthly satellite climatology of SST and Chi, and climatologic winds, while those in Table 3 are 
from monthly interannual satellite products and winds. The GoM and SS have the largest 
interannual variability in sea-air CO 2 flux. The flux in the SS is positive (source) in 2005 (+0.15 
Tg C yr' 1 ) and negative (weak sink) in 2006 (-0.02 Tg C yr' 1 ), while the largest flux (-1.55 Tg C 
yr' 1 ) occurred in 2007. These large differences in the SS annual fluxes are a result of large 
interannual changes in the spring drawdown of surface ocean pCOi (see Fig. 13). However, in 
the GoM the large differences in annual flux (+0.17 Tg C yr' 1 in 2004 and -0.19 Tg C yr' 1 in 
2007) are a result of wind speed variability as there are not significant interannual changes in the 
surface ocean pCOi seasonal cycle, as shown in Fig. 13. 

Averaged over the entire eight years, the MAB, GB+NS, and SAB are relatively the 
largest sinks of CO 2 to the atmosphere (-2.1, -1.0, and -0.9 Tg C yr' 1 , respectively), while the 
GOM is a small source (+0.01 Tg C yr' 1 ) and the SS a relatively small sink (-0.6 Tg C yr' 1 ), 
albeit with large changes from year to year. The east coast uptake (mean over the 8 years) is 4.6 
Tg C yr' 1 , which is at the upper end of the estimates from the binned field measurements with 
two different gas transfer parameterizations (4.0 to 4.3 Tg C yr' 1 ), and 3.6 to 4.0 Tg C yr' 1 from 
the algorithm using monthly climatology inputs (see Table 2). Table 3 shows that the lowest 
estimate occur in 2006 (3.4 Tg C yr' 1 ) and the highest in 2007 (5.4 Tg C yr' 1 ). 

The interannual variability in sea-air flux in all regions is mostly due to changes in the 
surface ocean pCOi, mainly in response to changes in solubility and biological drawdown due to 
variability in SST and phytoplankton production, respectively, and the wind-dependent gas 
exchange at the sea-air interface, accounted for by the gas transfer coefficient /ceeo (in cm hr' 1 ). 
From Table 1 we see that the algorithm pC02 sensitivity to the input parameters varies 
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significantly from region to region. In fact, the coefficients of many parameters change sign on a 
regional basis. So, in order to evaluate which parameters influenced the resulting estimates of 
sea-air flux the most, one needs to examine the yearly changes of these parameters and evaluate 
how much influence they have on the pCQi. As an example, there was a significant shift in the 
mean annual sea-air flux in the SS from 2005 to 2007 (Table 3 and Fig. 14). In 2005 the SS was 
a weak source of atmospheric CCF (+0.15 Tg C yr' 1 ), while in 2007 it shifted to a relatively 
strong CO 2 sink (-1.55 Tg C yr' 1 ). This shift was associated with lower SST (-0.8 °C), higher 
logio[C/7/] (+0.067), and higher k^o (+2.19 cm hr' 1 ) on average in 2007 compared to 2005. Using 
the coefficients for SS in Table 1, 8.77 ± 0.26 patm (°C)' 1 , -100.32 ± 4.66 patm (logioCC/?/])' 1 , we 
get the following changes in pC02 in 2007 compared to 2005: -7.1 ± 0.2 patm from SST and -6.7 
± 0.3 patm from Chi, for a total decrease in surface ocean pC02 of -13.8 ± 0.4 patm. Considering 
that this is a regionally and annually averaged value, this is a significant change in pCOj , which, 
combined with the increase in is the main reason leading to changes in sea-air flux. 

Time series (2003-2010) of annual mean sea-air CCA flux averaged for each of the five 
regions, each combined with annual means of SST, logio \Chl\, and k 66 o, are shown in Fig. 14. 

We show lo g 1 0 [ C/ 7 /] instead of absolute Chi concentration because the log-transformed Chi is the 
parameter used by the algorithm. Examination of each of these time series reveals some 
interesting interannual changes. The scale of variability for each variable changes from region to 
region and it is reflected by adopting different vertical axis ranges for each region. Interestingly, 
2006 is a year of transition for all regions north of Cape Hatteras (MAB, GB+NS, GoM, and SS). 
In 2006, the highest SST and Chi occur in the GoM and SS, followed by a decrease in SST 
reaching a minimum in 2007, which, combined with a peak in k 6 6o resulted in the largest uptake 
of CO 2 by the ocean in these two regions. As a result, there was a transition in the sea-air flux in 
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the SS from a very weak si nk in 2006 (-0.02 Tg C yr' 1 ) to a stronger sink in 2007 (-1.55 Tg C 
yr' 1 ). There was an increase of SST from 2007 to 2010 that contributed to a reduction in the 
ocean uptake. The sea-air flux interannual variability in the GB+NS, MAB, and SAB was also 
largely driven by changes in SST, with warmer years having reduced ocean uptake and colder 
years showing an increase in uptake. 

The annual mean time series of sea-air flux for each region (2003-2010), and the total for 
the entire east coast, are shown in Fig. 15. The GoM and SS regions were relatively stronger 
si nk s of CO 2 to the atmosphere in 2007 (-0.19 and -1.55 Tg C yr' 1 , respectively). The annual 
uptake of CO 2 ranged from -0.51 to -1.12 Tg C yr' 1 in the SAB with a mean of -0.89 ± 0.18 Tg C 
yr' 1 for 2003-2010. The equivalent values for the GB+NS were similar, with a range of -0.73 to 
-1.20 Tg C yr' 1 and an overall mean of -1.00 ± 0.18 Tg C yr' 1 . The MAB was the largest si nk 
with values ranging from -1.73 to -2.43 Tg C yr' 1 , and an overall mean of -2.12 ± 0.24 Tg C yr' 1 . 
The total sea-air flux (sum of all five regions) ranged from -3.4 to -5.4 Tg C yr' 1 , with the lowest 
uptake in 2006 and the highest in 2007. 

5.0 Summary and Future work 

We reconstructed a monthly climatology of surface ocean 77 CO 2 for the North American 
east coast continental shelf and developed regional algorithms to analyze the seasonal and 
interannual variability of surface ocean pCOi and sea-air CO 2 flux. A sensitivity analysis of 
parameters that influence the surface ocean pCOi showed that changes in DIC and SST are the 
main drivers for the pCOi seasonal cycle. Vertical mixing, mixing of low-salinity waters with 
shelf water, and biological drawdown are highly influential in the DIC variability. Much larger 
seasonal cycle amplitudes of DIC occur in regions north of Cape Hatteras than south of it. The 
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annual sea-air CO 2 flux for the entire East Coast derived from the algorithm ranges from -3.4 Tg 
C yr' 1 (2006) to -5.4 Tg C yr' 1 (2007) during the analyzed period (2003-2010). In general, 
estimates from the binned data and algorithm are in agreement with previous studies when the 
range of uncertainty and interannual variability are taken into account. 

Uncertainties in the estimates of sea-air flux can be reduced by filling the spatial and 
temporal gaps in the existing surface ocean /?COi inventory for the US East Coast. The 
limitations of spatial and temporal surface ocean pCC >2 data coverage present a challenge in 
validating algorithms and biogeochemical model pC02 and sea-air flux estimates. Improvements 
can only be obtained by continuous monitoring of pC 02 and other carbon cycle related variables 
in the near shore and shelf regions of the US East coast. As shown in Fig. 3, all regions have 
major spatial and temporal gaps in the data coverage. 

In this study, we used a multiple regression approach to convert regional satellite 
observed quantities (SST, and Chi) into pCOi . However, the relationship /?CC> 2 =/(SST, Chi, 
SSS, time) is empirical and does not represent a unique solution as pCOi depends on factors 
other than local SST and Chi, for instance. Surface waters with identical SST and Chi can 
possibly have different pC02 levels. However, there have been studies that apply the technique 
of neural networks for mapping in situ /AXE data in the open ocean [ Friedrich and Oschlies, 
2009; Lefevre et al., 2005; Telszewski et al., 2009]. The advantage of the neural network 
approach is that it can recognize and exploit relationships in the data which are not pre-defincd 
(as in regression techniques) and need to be expressible by an equation. This makes neural 
networks particularly suited to mapping relationships that are non-linear and empirical, provided 
sufficient data are available to ‘train’ the network. This technique looks promising for mapping 
the surface ocean /AO 2 in other coastal regions as well. 
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Hales et a!. [2012] presented a method for predicting coastal surface-water pC02 from 
remote-sensing data, based on self organizing maps (SOMs) and a nonlinear semi-empirical 
model of surface water carbonate chemistry, a method potentially applicable to the coastal 
regions in this study. The SOM approach was used to objectively map the sub-regions, while an 
entirely different approach was used to develop the pCC >2 algorithm within the SOM-defined 
sub-regions. The model used simple empirical relationships between carbonate chemistry (DIC 
and Aik) and satellite data (SST and Chi). Surface-water pCOi was calculated from the 
empirically-predicted DIC and Aik. This directly incorporated the inherent nonlinearities of the 
carbonate system, in a completely mechanistic manner. 
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Appendix A - Additional Sources of Surface Ocean pC0 2 not included in the SOCAT Data 

A. 1 South Atlantic Bight 

Underway surface ocean pC0 2 data from the SAB were collected by Dr. Wei-Jun Cai (a 
co-author in this study) and co-workers at the Department of Marine Sciences, University of 
Georgia. A total of 65,454 underway surface ocean pC0 2 records were processed for this study 
from six cruises along the SAB continental shelf: 5-16 January 2005, 19-30 March 2005, 27 
July to 5 August 2005, 7-17 October 2005, 16-21 December 2005, and 17-27 May 2006. The 
SOCAT data set includes the 2005 cruises but not those undertaken in 2006, which were added 
to our analysis to include all cruises. In all of the sampling cruises except for the one in 
December 2005, the research vessel transected the whole SAB from coastline to about 500-m 
water depth. The survey focused on 5 cross-shelf transects that are named E-, D-, C-, B-, and A- 
transect, respectively from north to south. In December 2005, the ship transected the whole SAB, 
but did not cover D- and B-transects and did not go beyond the 200 m isobaths due to limited 
ship time. Surface water and atmospheric xCC >2 were measured underway during all cruises. Sea 
surface temperature (SST) and salinity were recorded continuously with an onboard SeaBird 
flow through thermosalinograph. Sea level pressure was recorded using an onboard R.M. Young 
barometric pressure sensor. Surface water xCC >2 was measured using a LI-COR 7000 infrared 
gas analyzer coupled to a gas-water equilibrator. Details of the methodology and accuracy of 
instruments used are given in Jiang et al. [2008]. Fig. A-l shows the data distribution map. 

A.2 Gulf of Maine 

Underway surface ocean pC0 2 data from monthly cruises in the southern Gulf of Maine 
were obtained from the University of New Flampshire (UNIT) and integrated with the SOCAT 
data base. Underway data are measured continuously from pumped surface water for physical, 

33 


722 

723 

724 

725 

726 

727 

728 

729 

730 

731 

732 

733 

734 

735 

736 

737 

738 

739 

740 

741 

742 

743 

744 


chemical, biological and bio-optical properties. The data used in this study consisted of 309,665 
surface observations spanning the period of 2004-2010. These data originate from the UNH 
Coastal Ocean Observing Center's Coastal Carbon Group, which is an interdisciplinary research 
team within UNH-EOS engaged in efforts to observe and model how the Earth's pool of carbon 
moves between the land, ocean, and atmosphere with a particular focus on how this carbon 
cycling occurs in coastal regions, such as the Gulf of Maine. Dr. Joe Salisbury, a co-author in 
this study, is a member of the UNH Coastal Carbon Group. The methodology and 
instrumentation details are given in Vandermark et al. [2011]. The precision of the /CO? 
measurements was ±3 palm. Fig. A-2 shows the data distribution map. All underway cruise 
tracks are in the GoM, except for a single cruise track from Woods Hole to New York City. 

A. 3 Scotian Shelf 

Underway (UW) surface ocean pC0 2 data from transects across the Scotian Shelf, and 
high frequency pC0 2 , SST, SSS and calibrated fluorometer Chi data from the CARIOCA buoy 
were obtained from Dalhousie University [ Shadwick et al., 2010; Shadwick et al., 2011]. These 
data were used to evaluate the algorithm perfonnance on the Scotian Shelf. Hourly, autonomous 
observations of surface water pC0 2 (patm), chlorophyll-^ fluorescence (F C hi), and SST, were 
made using a CARIOCA buoy moored roughly 30 km offshore from Halifax, at 44.3° N and 
63.3°W, between April 2007 and June 2008. Hourly CARIOCA data were uploaded and 
transmitted daily via the ARGOS satellite system. The pC0 2 measurements were made by an 
automated spectrophotometric technique. A Sea-Bird (SBE 41) conductivity and temperature 
sensor was used to measure temperature (°C) and to determine salinity; chlorophyll-a 
fluorescence (pg l' 1 ) was determined by a WET Labs miniature fluorometer (WETstar). Non- 
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photochemical effects that are related to the intensity of the incoming solar radiation may 
decrease Fcm up to 80% during the day. This effect can be avoided by using night-time data 
which, to a large extent, are free of the effects of non-photochemical quenching, for fluorometer 
calibration. Night-time data were taken as a mean Fcm between 03:00 and 06:00 UTC (or 1 1 :00 
and 02:00 LT); data points were temporally interpolated to match discrete chlorophyll-a 
measurements (Chl-a in mg m' ) from monthly or twice monthly occupations at the mooring site. 

Chi -a concentration was determined fluorometrically in a Turner Designs fluorometer using the 

2 

acid ratio technique for seawater samples collected at 3, 5, or 10m depth. A linear regression (r 
= 0.76, N=29 ,p < 0.001) was used to determine the relationship between the Fchi and Chi -a, and 
applied to the CARIOCA fluorescence-derived Chi -a time-series (Chi/ in mg m' ). Shadwick et 
al. [2010] perfonned a validation of satellite monthly chlorophyll data by regressing it against 
the (night-time calibrated), monthly mean, CARIOCA Chi/.- time series, (r =0.68, N=1 4, p 
<0.002). 

Measurements of pC0 2 UW were made by a continuous flow equilibration system in: 
October 2006, April, August, and October 2007, and April and October, 2008 on board the 
CCGS Hudson. The UW measurements (see distribution map in Fig. 7a) were obtained on 
monitoring cruises on the Scotian Shelf (see Shadwick et al., 2011 for details of the field 
program). Measurements of pC0 2 UW were made by a non-dispersive, infrared spectrometer 
(LiCor, LI-7000). The system was located in the aft-laboratory of the ship and the intake depth 
was approximately 3m below the water surface. Measurements were made every minute and 
used to compute hourly averages. The system was calibrated daily with both a CCfi-free 
reference gas (N 2 ) and a CO 2 calibration gas (328.99 ppm) provided by the US National Oceanic 
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and Atmospheric Administration (NOAA). The data were corrected to in-situ water temperature 
and to 100% humidity and had an associated uncertainty of less than lpatm. 
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Table 1. Coefficients and statistical data for pC0 2 algorithm (Equation 1). The multiple regression coefficients and 
their corresponding standard errors were obtained using the MatLab function ‘regstats’ with t statistics. 


Means and Coefficients 

SAB 

MAB 

GB+NS 

GoM 

SS 

T 0 (°C) 

23.21 

15.27 

11.27 

10.29 

7.34 

So (psu) 

35.38 

31.64 

32.19 

31.41 

30.58 

Y (days) 

123 

218 

359 

343 

27 

Chl 0 (mg m' 3 ) 

1.09 

1.54 

1.62 

2.94 

1.24 

a (patm) 

378.69 

±1.76 

360.07 

±1.40 

370.66 

±1.84 

373.06 

±1.38 

351.43 

±0.90 

b (^atm) 

24.00 

±2.05 

7.03 

±4.82 

37.05 

±2.63 

39.43 

±1.68 

69.31 

±2.39 

c (patm “C' 1 ) 

12.23 

±0.36 

5.20 

±0.47 

6.88 

±0.40 

1.65 

±0.24 

8.77 

±0.26 

d (jaatm psu’ 1 ) 

-22.49 

±1.71 

1.11 

±0.61 

-10.95 

±2.33 

-1.34 

±0.83 

1.44 

±0.86 

e (|iatm/log 10 (Cb/)) 

30.25 

±5.87 

-14.99 

±5.51 

10.05 

±7.67 

-20.65 

±3.83 

-100.32 

±4.66 

2 

r 

0.82 

0.55 

0.60 

0.42 

0.74 

RMSE (fjatm) 

26.7 

36.9 

32.2 

34.6 

22.4 

N 

356 

997 

356 

847 

684 
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Table 2. Sea-air CO 2 flux for reference year 2004 from binned data, algorithm for year 2004, and previous studies 
(literature). Uncertainties (e) for the estimates from the data were obtained as e = STD/VN, where STD is the 
standard deviation and N the number of data points. Uncertainties for the sea-air C0 2 flux estimates from the 
algorithm were based on the standard deviation of all monthly estimates for the period 2003-2010. Both specific 
(mol C0 2 m" 2 yr' 1 ) and total (Tg C yr' 1 ) sea-air fluxes are shown for each region and total for the whole coast. Two 
gas transfer coefficients were used, the polynomial equation of Wanninkhof et al. [2009] (fcg 60 ) and the quadratic 
dependence version of Ho et al. [2011] (fc| 6 o) adjusted for steady winds using the nonlinearity coefficients C 2 and 


C 3 . 


Region 

Area 
10 10 m 2 

Data 

mol C0 2 m" 2 yr" 1 
Tg C yr- 1 

Algorithm 
mol C0 2 m" 2 yr" 1 
Tg C yr- 1 

Literature 
mol C0 2 m" 2 yr" 1 
Tg C yr- 1 



^660 

^660 

^660 

^660 


SS 

12.82 

-1.10 + 0.25 

-1.21 + 0.27 

-0.39 + 0.34 

-0.42 + 0.36 

+1.42 + 0.28 (d) 



-1.69 + 0.39 

-1.87 + 0.42 

-0.56 + 0.50 

-0.60 + 0.53 

+2.19 + 0.43 

GoM 

12.77 

+0.11 + 0.21 

+0.04 ± 0.22 

+0.01 + 0.08 

+0.01+0.08 

+0.38 ± 0.26 (c) 



+0.17 + 0.32 

+0.06 + 0.34 

+0.02 + 0.12 

+0.02 + 0.12 

+0.58 ± 0.40 

GB+NS 

5.83 

-0.65 + 0.20 

-0.71 + 0.22 

-1.27 + 0.23 

-1.37 + 0.24 

- 



-0.46 + 0.14 

-0.50 + 0.15 

-0.79 + 0.16 

-0.86 + 0.16 

- 

MAB 

9.31 

-0.95 + 0.24 

-1.07 + 0.27 

-1.58 + 0.19 

-1.78 + 0.19 

-1.1 ±0.7 



-1.06 + 0.27 

-1.12 + 0.30 

-1.63 + 0.21 

-1.83 + 0.22 

-1.0 ±0.6 (a) 

SAB 

10.20 

-0.79 + 0.26 

-0.68 + 0.24 

-0.61 + 0.17 

-0.67+0.16 

-0.48 + 0.21(6) 



-0.97 + 0.31 

-0.83 + 0.29 

-0.67 + 0.20 

-0.74 + 0.20 

-0.59 + 0.26 

Total 

50.63 

-4.01 + 0.30 

-4.26 + 0.31 

-3.63 + 0.24 

-4.01 + 0.25 

- 


(a) DeGrandpre [2002]; ( b ) Jiang et al. [2008]; (c) Vandemark et al. [2011] is 5-year mean (2004-2208) 
but ranging from +0.71 (2005) to -0.1 1 (2007) mol m" 2 yr" 1 ; ( d) Shadwick et al. [2011]. Values for (b), 
(c), and (d) were converted from specific to total flux, or mol C0 2 m" 2 yr" 1 to Tg C yr" 1 (xl2xareaxl0 12 ). 
k\ b 0 = 3 + 0.1 U 10 + 0.064C 2 U? o + 0.011 C 3 Ul 0 and k 2 660 = 0 ,262C 2 U? 0 
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Table 3. Sea-air C0 2 flux derived from the regional algorithms for 2003-2010. The flux is given in two different 
units for each year (mol C0 2 m' 2 yr _1 /Tg C yr" 1 ), and in Tg C yr' 1 for the overall 8-year mean and whole coast sum. 
The flux was calculated using the gas transfer equation of Ho et al. [2011]. 


Year 

SAB 

MAB 

GoM 

GB+NS 

SS 

Sum 

2003 

-0.78/-0.90 

-2.18/-2.43 

+0.002/+0.009 

-1.72/-1.20 

-0.33/-0.55 

-5.07 

2004 

-0.75/-0.88 

-2.08/-2.31 

+0.107/+0.166 

-1.72/-1.20 

-0.27/-0.39 

-4.61 

2005 

-0.95/-1.12 

-1.92/-2.13 

+0.068/+0.108 

-1.49/-1.04 

+0.18/+0.15 

-4.03 

2006 

-0.74/-0.88 

-1.56/-1.73 

-0.052/-0.074 

-1.05/-0.73 

-0.01/-0.02 

-3.43 

2007 

-0.43/-0.51 

-1.76/-1.95 

-0.129/-0.191 

-1.71/-1.20 

-1.01/-1.55 

-5.40 

2008 

-0.78/-0.93 

-1.72/-1.91 

-0.045/-0.062 

-1.21/-0.85 

-0.55/-0.77 

-4.52 

2009 

-0.66/-0.76 

-1.90/-2.1 1 

-0.024/-0.028 

-1.32/-0.92 

-0.72/-1.14 

-4.96 

2010 

-0.91/-1.08 

-2.16/-2.41 

+0.079/+0.126 

-1.21/-0.85 

-0.18/-0.40 

-4.62 

Mean 

-0.89±0.18 

-2.12±0.24 

+0.007±0.1 12 

-1.00±0.18 

-0.58±0.52 

-4.58 
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Fig. 1. Regional domains for analysis adapted from Hofmann et al. [2008] and Hofmann et al [201 1], The white 
circles show the locations of the NDBC buoys within each regional domain. The white star shows the location of the 
Sable Island meteorological station and the white square the location of the Carioca buoy. 
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Fig. 2. Seasonal climatology maps of SST, SSS, and Chi. Upper row: SST composites from MODIS Aqua; middle 
row: SSS composites from World Ocean Data 2009; bottom row: Chi composites from MODIS Aqua. Refer to the 
methods section (3.0) for details. The MODIS SST and Chi seasonal climatologies are based on the period 2002- 
2011. The seasons are defined as Dec-Jan-Feb (DJF), Mar-Apr-May (MAM), Jun-Jul-Aug (JJA), and Sep-Oct-Nov 
(SON). 


2 



300 350 400 

pC0 2 (f^atm) 


450 500 0 1 2 3 4 5 6 7 8 9 10 11 12 

Number of Months 


Fig. 3. Color-coded SOCAT surface ocean pC0 2 cruise tracks (a) and corresponding coastal binned data (b) with 
associated color-coded temporal coverage in months. The highest temporal coverage corresponds to the most 
travelled routes (in orange to red), i.e., most frequent destination ports (Boston, New York, Norfolk, Miami) used by 
the Volunteering Observing Ships. The SOCAT data set also includes transects occupied by research vessels. The 
SS, GoM, GB+NS, MAB and SAB regional boundaries are overlaid as black lines. 
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Fig. 4. Plot of goodness-of-fit statistics for all regional MLRs with incremental addition of corresponding proxy 
parameters. The x-axis shows the RMSE normalized by the maximum attained value among all MLRs, while the 
axis shows (1 -r 2 ). Thus a perfect match between data and MR values would be centered at the origin (0, 0). 
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Fig. 5. From top to bottom: scatter plots (left column) of observed (SOCAT) vs. algorithm (Equation 1) pC0 2 
(patm) for the five regions (black dots all months, green squares monthly ensemble averages). The right column 
shows the mean seasonal plots of the ensemble averages for the equivalent regions. There are no data available for 
the MAB and GB+NS for January. Only data bins with more than six months of coverage were used. 
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Fig. 6. (a) Scatter plot of algorithm versus observed surface ocean pC0 2 based 
on observed values not used in the algorithm development (bins with temporal 
coverage less than 6 months). The r 2 , RMSE, and mean absolute percent 
difference (MAPD) are shown in the legend, (b) Histogram of residuals (observed 
minus algorithm). The red dashed vertical lines represent the standard deviation (±ct) 
of the observed pC0 2 . 
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Fig. 7. (a) Map showing the seven 2°x2° boxes covering the entire Scotian Shelf (SS) region adapted from Shadwick 
et al. [2010]. The contour line is the 200 m isobath. The algorithm and in situ (SOCAT (not shown) and UW 
observations from Dalhousie University cruises) mean surface ocean pC0 2 were obtained for each of the seven 
boxes for evaluation purposes. The scatter plot of algorithm versus observed pC0 2 for all seven boxes is shown in 
(b) with corresponding statistics in the legend. 
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Fig. 8a. Time series of algorithm mean surface ocean pC0 2 (black lines) for boxes 1 through 4 shown in Fig. 7a. The 
corresponding SOCAT (red dots) and Dalhousie UW (blue dots) data are shown for comparison. The blue lines are 
the atmospheric pC0 2 . See Fig. 7b for statistical evaluation. 
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Fig. 8b. Time series of algorithm mean surface ocean pC0 2 (black lines) for boxes 5, 6 and 7 shown in Fig. 7a. The 
corresponding SOCAT (red dots) and Dalhousie UW (blue dots) data are shown for comparison. The blue lines are 
the atmospheric pC0 2 . See Fig. 7b for statistical evaluation. 
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Fig. 9. Time series of high frequency (hourly) surface ocean pCOi measured (blue crosses) at the Carioca buoy from 
2007 to 2010, and corresponding algorithm prediction (red crosses) using hourly values of SST, SSS, and calibrated 
fluorometer Chi as inputs (top panel). The scatter plot of observed vs. algorithm pC0 2 is shown in the bottom panel. 
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Fig. 10. Seasonal maps of algorithm pC0 2 , salinity-derived alkalinity from Cai et al. [2010] equations, and DIC 
derived from alkalinity and algorithm pC0 2 - The seasons are defined as Dec-Jan-Feb (DJF), Mar-Apr-May (MAM), 
Jun-Jul-Aug (JJA), and Sep-Oct-Nov (SON). 
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Fig. 11. Sensitivity ofpCO? seasonal cycle to most influential parameters. Alkalinity was derived using SSS from 
monthly WOA 2009 salinity data (D. Tomaso personal communication, 2012), spatially interpolated using Kriging, 
and Cai el al. [2010] equations. DIC was derived from algorithm pC0 2 , alkalinity, WOA SSS, and MODIS SST. 
Refer to text for methodology to derive parameter sensitivity. 
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Fig. 12. Regionally averaged seasonal DIC (black lines and circles) derived from TA (SSS) [Cai et al., 2010], SST 
from MOD1S, monthly SSS from WOA 2009 ( D . Tomaso personal communication, 2012) spatially interpolated 
using Kriging, and algorithm pC0 2 - The seasonal mixed layer depth (MLD) is superposed for each region (red lines 
and circles). The red dashed lines represent the mean bottom depth for each region and the thin black lines are the 
annual mean DIC for each region, with the GoM and SAB having the highest values (2022 pmol kg" 1 ) and the MAB 
the lowest (1968 pmol kg' 1 ). 
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Fig. 13. Left panel: Monthly surface ocean pC0 2 derived from algorithm (black lines) and atmospheric pCXh from 
Griffon, NC located at 35.53°N and 77.38°W (superposed blue lines). Right panel: Sea-air CCL flux derived from 
ApC0 2 , CCMP winds, and Ho et al. [2011] gas transfer parameterization. 
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Fig. 14. Mean annual sea-air C0 2 flux (red lines, Tg C yr" 1 ) combined with SST (°C), loglO[CA/] (blue lines) and 
k 660 (cm hr' 1 , blue lines) for all 5 regions 
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Fig. 15. Time series of algorithm annual sea-air C0 2 flux for all five individual regions and for the entire east coast. 


16 



17 


600 


550 
500 
450 
400 
350 
300 
250 
200 
150 
100 

Fig. A-2. Map showing the underway pC0 2 tracks in the GoM and a single cruise track from 
Woods Hole to New York City. 
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